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Abstract 

We bosonize a Fermi liquid in any number of dimensions in the limit of long 
wavelengths. From the bosons we construct a set of coherent states which 
are related with the displacement of the Fermi surface due to particle-hole 
excitations. We show that an interacting hamiltonian in terms of the original 
fermions is quadratic in the bosons. We obtain a path integral representation 
for the generating functional which in real time, in the semiclassical limit, gives 
the Landau equation for sound waves and in the imaginary time gives us the 
correct form of the specific heat for a Fermi liquid even with the corrections 
due to the interactions between the fermions. We also discuss the similarities 
between our results and the physics of quantum crystals. 
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I. INTRODUCTION 



Attempts to describe fermionic systems by bosons date to the early days of second 
quantization. In the early 50's Tomonaga 0], generalizing earlier work by Bloch M on 
sound waves in dense Fermi systems, gave an explicit construction of the Bloch waves for 
systems in one dimension. Three years after the seminal work of Tomonaga, Bohm and 
Pines || showed that there was a natural connection between bosonization and the Random 
Phase Approximation (RPA). Subsequently many authors ||] derived an explicit Fermi- 
Bose transmutation in one-dimensional systems. These works uncovered deep connections 
in relativistic field theories (both fermionic and bosonic) and with condensed matter systems. 

The success of the bosonization approach in one dimension is related to phase space con- 
siderations. Even for non-interacting fermions, two excitations with arbitrarily low energies 
moving in the same direction move at the same speed (the Fermi velocity) and, hence, are 
almost a bound state. Consequently, even the weakest interactions can induce dramatic 
changes in the nature of the low-lying states. These phase space effects become manifest 
even in perturbation theory and result in the presence of both marginal and marginally 
relevant operators. Bosonization yields a simple and very general description of the possi- 
ble fixed point structure of strongly correlated systems: marginally relevant operators are 
responsible for the gap generating instabilities while strictly marginal operators give rise to 
anomalous dimensions. Quite generally, the fixed points are determined by the marginal 
operators alone. In turn, the anomalous dimensions produce non-Fermi liquid behavior. 
Bosonization in one spacial dimension is the statement that the charge and current densities 
for states restricted to the vicinity of the two Fermi points, p(x) and j(x), of electrons obey 
the equal-time algebra (Kac-Moody) 

\p{x)J(y)] = -~6(x-y) (1.1) 

This algebra implies that there exists of a free bosonic field (f)(x) and its canonically conjugate 
momentum IT(a;) which obey the canonical equal-time commutation relation 
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[00*0, n(y)] =%5{x-y) 



(1.2) 



provided that one makes the identifications 



p(x) 



1 d 



(f)(x), j(x) = -7=II(a;) 



(1.3) 



More complicated operators, such as various order parameters associated with 2k f excita- 
tions, can also be identified with suitably chosen bosonic operators ||. Typically, these 
operators involve non-local or non-linear (or both) expressions in the bosonic field. In any 
event, the basic building blocks are the local density and current operators which describe 
electron-hole pairs at small momentum transfers. 

In dimensions higher than one, phase space considerations change the physics of the low- 
lying states. The only marginal operators left (i. e. , operators whose coupling constants 
(Landau parameters) do not change as the energy scale is lowered) are responsible only for 
changing the shape of the Fermi surface. However the low temperature and low frequency 
behavior of physical quantities is insensitive to the presence of such marginal operators. In 
a sense, these operators are redundant. This observation is at the root of the stability of the 
Landau theory of the Fermi liquid 0,0. It is, thus, hardly surprising that very few attempts 
have been made to generalize the bosonization approach to dimensions higher than one. In 
addition, phase space considerations also tell us that even for small momentum transfers, an 
electron-hole pair can decay into its electron and hole constituents. Electron-hole collective 
modes do exist ( and are described by the Random Phase Approximation) but so does the 
electron hole continuum. Hence, unlike what happens in one space dimension, kinematics 
alone does not require the generic existence of a bosonic bound state. 

The first serious attempt at bosonization in higher dimensions was carried out by Luther 
H who constructed a generalized bosonization formula in terms of the fluctuations of the 
Fermi sea along radial directions in momentum space. However this approach was not 
pursued. Interest in the construction of bosonized versions of Fermi liquids has been re- 
vived recently in the context of strongly correlated systems due to the possibility of 
non-conventional ground states which are not Fermi liquid like. Since in one dimension 
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the bosonization approach is a powerful non-perturbative tool, the expectation is that if a 
suitably generalized analog does exist in higher dimensions, it may throw light on the non- 
perturbative aspects of Fermi liquids and, hopefully, it may help to define non- Fermi liquids. 
It is well known since the early days of the Landau theory that particle-hole excitations have 
bosonic character e. g. , the sound waves (zero sound collective modes) of the Fermi surface 
of neutral liquids or plasmons in charged Fermi liquids || . Inspired by these identifications, 
and drawing from Luther's work, Haldane |10| has derived recently a bosonic algebra for 



density fluctuations of a Fermi liquid in the form of a generalized Kac-Moody algebra which 
generalizes eq (1.1). Haldane's construction has been examined recently by several authors 

H0- 



Our approach has the same starting point as Haldane's. We also begin by constructing 
a bosonic algebra with creation and annihilation operators which act over a reference state 
describing a filled Fermi sea. This algebra is valid in a restricted Hilbert space of states 
generated by small deformation of the filled Fermi sea. Using these operators we construct a 
set of coherent states which can be regarded as deformation of the Fermi surface. More pre- 
cisely, we show that the eigenstate of the annihilation operators of the algebra, namely the 
coherent states, are related with the displacement of the Fermi surface in some direction. 
These displacements are a coherent superposition of particle-hole excitations close to the 
Fermi surface. This interpretation leads to the picture of the Fermi surface as a real dynam- 
ical object. We have two strong arguments in support of this picture. The first argument 
is based on the fact that, using these coherent states, we reproduce the phenomenological 
theory of Pomeranchuk || on the stability of a Fermi liquid and with its modern reinter- 
pretation (by Shankar j7|) in terms of the renormalization group method. Pomeranchuk's 
approach is based on the work of Landau on many-body systems. The second argument fol- 
lows from the fact (which we elaborate further below) that it is possible to define a surface 
tension associated with the Fermi surface itself. 

Armed with these coherent states, we construct a coherent state path integral which can 
be viewed as a sum over the histories of the shape of the Fermi surface, a bosonic shape field. 



We show that there exists a class of simple interacting hamiltonians (in terms of the original 
fermions of the theory), which are naturally related to the Landau theory of the Fermi liquid, 
which in the coherent state bosonization leads to a simple quadratic hamiltonian in terms 
of the bosons. These hamiltonians are the natural generalization of the Luttinger/Thirring 
models of one dimensional systems. According to the renormalization group analysis, these 
hamiltonians represent the stable fixed points of Fermi liquids. 

With the path integral we can study the thermodynamics and dynamics of interacting 
electronic systems in a very simple way. In order to prove the consistency of our method we 
calculate the low temperature specific heat and obtain the same results of the Fermi liquid 
theory (although we have only real bosons in the problem). We can go even further and 
calculate the corrections for the specific heat due to interactions and we show that our results 
also agree with the well known results of the Fermi liquid theory. This calculation reveals 
the new features of the bosonic field involved in the problem. They are not an usual free 
bosonic field but a topologically constrained excitation. The constraint comes from the need 
to sum over the bosonic excitations tangent to the Fermi surface and do not contribute to 
the energetics. This is most important since free non-relativistic bosons in space dimensions 
larger than one do not have the linear specific heat which is characteristic of Fermi liquids. 

In this paper we discuss the bosonization approach to the fixed point hamiltonians. 
These hamiltonians contain operators which are at most marginal. We deliberately leave 
out a number of relevant operators which lead to the well known instabilities of the Fermi 
liquid: superconductivity, magnetism, etc. We will discuss these operators elsewhere. 

We also show that the Landau theory is the semiclassical approximation for the bosons 
which is exact for the fixed point hamiltonians since the hamiltonian is quadratic in bosons. 
This result means that the bosonic field is nothing but the field of sound waves which prop- 
agates on the Fermi surface. We applied these methods to study the physics of sound waves 
in two dimensions. We obtain a general equation for the sound modes and, in particular, 
we study the zero sound. In order to study the characteristics of these sound modes we go 
to the limit of large Fermi momentum and concentrate in the forward scattering direction. 
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We introduce a new scale in the problem, the range of interaction. We find a rather com- 
plex spectrum of sound modes which describe increasingly wrinkled Fermi surfaces. In the 
case of long range interactions, the shape of the Fermi surface can become unstable. The 
existence of strong local quantum fluctuations of the shape of the Fermi surface also sug- 
gests an analogy with the quantum mechanical fluctuations of the shape of a crystal at zero 
temperature. Thus it is natural to ask if it is possible that, for long range interactions, the 
Fermi surface could undergo a roughening transition. In the case of of the shape of a three 
dimensional crystal it is known [14] that quantum fluctuations generally make the surface 
smooth, not rough. However, in the case of a "planar crystal" , which has a one dimensional 
surface, it is possible to have a quantum mechanical roughening transition. It turns out, 
however, that by an explicit calculation of the correlation function between different pieces 
of the Fermi surface we can show that quantum fluctuations of the fixed point hamiltonian 
wash out this interesting possibility. The reason is that the dynamics described by the fixed 
point hamiltonians enforce Luttinger's theorem as a local condition. It may be possible to 
find other hamiltonians which enforce Luttinger's theorem only as a global constraint. The 
shape of the Fermi surface of the ground states of these hamiltonians may exhibit a quantum 
roughening transition in two space dimensions. 

The paper is organized as follows. In sections II and III we present two physical argu- 
ments which suggest that the Fermi surface should be regarded as a quantum mechanical 
object. In section II we review the stability of a Fermi liquid and specialize the discussion to 
the case of two dimensions . In section III we refine the analogy with the theory of surfaces 
by defining the surface tension of the Fermi surface. In section IV we describe the coherent 
state bosonization construction for an interacting fermionic system. In section V we show 
how a fermionic hamiltonian can be rewritten in terms of the bosons and in section VI we 
study the Fermi liquid properties, the thermodynamics and classical dynamics. Section VII 
contains our conclusions. 
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II. THE STABILITY OF THE FERMI LIQUID IN TWO DIMENSIONS 



The route for the stability of a Fermi liquid was established by Pomeranchuk || who 
studied the effect of the change of the Fermi surface in the free energy of the system using 
the Landau expansion for the change in the total energy, 

AE = J2 e^rv + i J2 ftfSngSnp, (2.1) 
p p,p' 

where V is the volume of the system, Snp is the deviation of the occupation number with 
p is at the Fermi surface, e~ is the bare dispersion relation and f^p is the quasi-particle 
interaction. 

In this section we essentially follow the work of Pomeranchuk |),T3] but we restrict 
ourselves to two dimensions since in three dimensions the results are well known and two- 
dimensional Fermi liquids are of direct physical interest. Consider a change in the Fermi 
momentum given by 

Pf(0) = Pf + S PF (6) (2.2) 

where 9 is the angle which parametrizes the position of the Fermi surface and p F is the 
original Fermi momentum. 

The new occupation number is simply given by 

ry= 9(^(0) -p), (2.3) 

substituting (2.2) in (2.3) and expanding up to second order in the deviations we find 

dnp = 5p F {6)5(p F -p) — . (2.4) 

Now, using (2.1) and (2.4) it is easy to show that the total change in the free energy, 
F = E — fiN (/x is the chemical potential and iV is the total number of particles ) is given 
by 

AF PF v F f^dO 2 p 2 F [**d9 f^dd' -j 

V Zn Jo Itx Its 1 Jo 2vr jo Ztt 
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where p and p' are at the Fermi surface and v p is the Fermi velocity. 

We expand the interaction and the displacement of the Fermi surface in Fourier compo- 
nents (in two dimensions we can parametrize the points of the Fermi surface by one angle, 
9 ) as, 

+00 

f(p,p') = J2fm™s(m(9-9>)) (2.6) 

m=0 

and 

+00 

Sp F (9) = u m cos(m9), (2.7) 

771=0 

substituting these expressions in (2.5) and using the orthogonality between the Fourier 
components one finds, 

where N(0) = - E£ - is the density of states at the Fermi surface. 
Observe that the Fermi liquid is stable if AF > or 

N(Q)f n > -1 (2.9) 

for all n. This result implies that the Fermi liquid is a local minimum in the configuration 
space of the many-body system. Local stability, in terms of renormalization group, means 
that the effective hamiltonian is at a fixed point [0. 

We notice that it is not a coincidence that the free energy (2.8) has the same form as 
the one for a drumhead. It means that the Fermi surface is a dynamical object with elastic 
properties. Furthermore, we see that the interaction term gives the same order of contribu- 
tion as the free term of the energy. In terms of the renormalization group interpretation of 
the stability of a Fermi liquid J7[ it means that the interaction is a marginal operator while 
next order expansion in the functional (2.1) will give only irrelevant operators which do not 
contribute to the low energy physics of the problem. 

Anderson |1 has argued that a singular interaction between the fermions could give rise 



to a new features which could explain the anomalous behavior of cuprates. Haldane |L0f has 



constructed an argument which is mainly based on the form of expression (2.8) which we 
reproduce here due to its simplicity and elegance. Suppose that instead of an overall density 
of states, N(0), we introduce a local density of states at each point of the Fermi surface. In 
order to do so we have to define a cut-off on the Fermi surface, A, which is much smaller 
than the Fermi momentum but yet large enough to count a macroscopic number of states. 
By naive dimensional analysis we conclude that this local density of states must scale like 
A d ~ 1 Vp 1 . Therefore, from (2.8), the quantity F m = A d ~ l Vp 1 f m is the dimensionless coupling 
constant. Suppose that the interaction, f m , and the Fermi velocity are well behaved (they 
do not diverge) in the limit where A — > (the scaling limit). It is easy to see that the F m 
goes to zero for d > 1 and therefore only the non- interacting term in present. This result 
implies that the interaction is marginally irrelevant. However, for d — 1, F m will be always 
finite, that is, the interaction is marginally relevant in one dimension. Therefore, even if the 
interaction is not singular in the scaling limit it has profound consequences in one dimension 
while it "scales away" in higher dimensions. However, we have the interesting possibility 
that the interaction term f m or the Fermi velocity diverge in the scaling limit. If this is so, 
the dimensionless coupling constant can diverge and therefore the operator becomes relevant 
with respect to the non-interacting term. However, this singular behavior is never present 
in a Fermi liquid since the renormalization group flow is dominated by the non-interacting 
fixed point. 

The approach in this section is essentially phenomelogical because we have postulated 
(2.1) in order to get the stability condition (2.9). In the next section we study the elastic 
properties of the Fermi surface from the microscopic point of view. This new point o view 
will bring a different perspective on the same subject, namely, the Fermi surface as a real 
quantum object. 
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III. THE TENSION OF THE FERMI SURFACE 



In this section we explore in more detail the picture of the Fermi surface as a quantum 
object. In particular we will establish that, for systems which obey the hypothesis of a Fermi 
liquid, it is possible to define a surface tension. Here we follow the work of Luttinger [13| 



on the properties of many-body systems. It can be shown that the total energy of a spinless 
interacting fermionic system can be written in the form [IB] of an integral over the whole 
momentum space as, 

£ = (3.1) 

k 



f+°° dw [w + e°r) 
E k = I ^ — n — G(k,w) (3.2) 



where 



-oo 2ni 2 

where i] — > + , is the non-interacting dispersion relation, G(k, w) is the interacting Green's 
function in momentum space. 



Equation (3.2) can be rewritten in terms of the spectral function of the system as [|17 



r+n r\ w (w + e? ) 
&i= I —^r^A(k\w) (3.3) 



oo 7T 



where the spectral function A(k, w) is given by 



A(k, w) = E f' w) 2 _ (3.4) 



w 



k 

where Si?(fc, w) and T,j(k, w) are the real and imaginary part of the self energy and fi is the 
chemical potential of the problem. 

The spectrum of the interacting system (the poles of the Green's function), e^, is given 
by the solution of the self consistent equation 

e E - 4 - £ a (£, ej) = (3.5) 

and, in particular, the Fermi surface is defined by the set of vectors {pf} such that 
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tp F = V- (3.6) 

It is clear from the above definitions that we expect that any singularity in this problem 
should appear near the Fermi surface. We, therefore, split the integral in (3.3) in two pieces, 
from — oo to \i — £ (free of singular terms) and from fi — £ to fi (the singular contribution) 
where ( is an small energy scale which permit us to look only to the singular part of Er. 

Moreover, consider the vectors k close to the Fermi surface, 

k = p F + 8k (3.7) 

where the displacement 8k is perpendicular to the Fermi surface, that is, 

8k.Ve^ F = 5k | Ve p - F | (3.8) 

and 8k is positive if the vector points outside of the Fermi surface and negative it points 
inside the Fermi surface. 

Since we are interested in the integral close to the Fermi surface we can expand the 
integrand using the above formulae. In particular the denominator in (3.4) can be rewritten 

as 

w-4- E R (k,w) = ZZ\ W - Cjf ) (3.9) 

where 



is the quasi-particle residue (T7 
We can also write 



e % = fj, + 8k | Ve p - F | (3.11) 



where we used (3.6), (3.7) and (3.8). 

Substituting these approximations in the singular part of the integral (3.3) and changing 
variables with respect to the chemical potential we get 
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* Jo tt 2 (zr 2 (w + 5k\ Ve p - F |) 2 + E 7 (fc, /x - w) 2 " 



For a Fermi Liquid we use the fact that the imaginary part of the self energy can be 
written as [|13|]- 



T>i(k,w) = C^w — fi) 2 sgn(fi — w) (3.13) 

where Cg depends only on k and sgn(x) = 1(— 1) if x > 0(x < 0). 

Substituting (3.13) in (3.12), changing the variables of integration from w to x = ^ Sk ^ tp F^ 
we found, 

(3.14) 

Observe now that if we let 5k — > we can use the well known expression 

Lim m ^ ^l^lir^L i V7 12Y = * 1 Zp > 1 (5(;E + s ^(^)) " S ( x )) ■ 

\Z^ x z (x + sgn{ok)Y + Cidk/ \ vtp F \ z J 

(3.15) 

Substituting this result in (3.14) we find 

EI=\Z fF \ Vpf 2 ^ e(-Sk) (3.16) 

where 6(x) = 1(0) if x > 0(x < 0). 

From (3.16) we conclude therefore that there is a discontinuity at the Fermi surface which 
is given by 

AE(p F ) =\ Z fF I Vpf 2 J (3.17) 

Therefore there is no singularity when the quasi-particle residue vanishes and we have 
non- Fermi liquid behavior. In the non-interacting system is easy to see that this discontinuity 
is simply given by \i which is the energy needed to put another electron in the system. 
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We can also show (see appendix A) that this result holds true if the imaginary part of the 
self energy goes like (w — with v greater than zero. However, if v is smaller than zero 

we can prove that the discontinuity in the E(k) vanishes as 5k^ l ~ v ^ 2 as 6k — > 0. Therefore, 



v = is a critical point (this is the case of a marginal Fermi liquid phenomenology [T8| ) . 

We can interpret this discontinuity in E(k) as due to a surface tension in momentum 
space since E{k) is the contribution of the mode k to the total energy, as we see by (3.1). 
The surface tension in three dimensions can be defined by the difference of energy density 
across the Fermi surface as follows 

aF [M (3 ' 18) 

where R F and R 2 F are the principal radii of curvature at pp and Vp is the volume of the 
Fermi sea. 

In two dimensions the analogous of (3.18) is 

where Ap is the area of the Fermi sea. 

Comparing (3.18) or (3.19) with (3.17) we conclude that the surface tension is propor- 
tional to the quasiparticle residue. Therefore, for a non-Fermi liquid behavior the surface 
tension vanishes with the quasiparticle residue. From this observation we can conclude that 
a transition from Fermi liquid to non-Fermi liquid behavior can be viewed as a phase tran- 
sition in momentum space for the case where the quasiparticle residue can be tuned to zero 
adiabaticlly. 

IV. BOSONIZATION AND COHERENT STATES 

Our starting point to approach the bosonization resembles the microscopic approaches for 
the foundations of Fermi liquid theory ]I5| . However, instead of working with the dynamics 
of the response functions, we will work directly with the properties of the operators which 



13 



generate the physicals spectrum, in a restricted Hilbert space of states of the filled Fermi 
sea. This is the standard procedure of bosonization in one dimension and in a recent work of 



Haldane | l0|| . From this perspective, the algebra obeyed by the operators is not a property 
of the operators themselves but a property of the states. 

For simplicity we will consider a system of interacting spinless fermions. Generalizations 
to systems with an internal symmetry, such as spin SU(2), can be done with some minor but 
important modifications. The density of fermions at some point f at some time t is given by 



p(r, t) = ^(r, t)m t) = £ 4-f (%+?(*) C 4 - 1 ) 

k,q 



where ct and are the creation and annihilation operator of an electron at some momentum 
k which obey the fermionic algebra, { c £' c fc'} = where {...} is the anticommutator 



and all other anticommutation relations are zero. In the Fermi liquid theory [15] the 
operator which appears in the r.h.s. of (4.1) determines the behavior of the system. We will 
concentrate on this operator which we will denote by 

%(M) = 4-f W^+fW- ( 4 - 2 ) 

In particular, rio(k) is the number operator in momentum space. As in the standard ap- 
proaches for the Fermi liquid theory we will concentrate our interest in regions close to 
the Fermi surface (k ~ pp) and consider long wavelength fluctuations around these regions 
(7-0). 

The equal time commutation relation between the operators defined in (2) is easily 
obtained. In the long wavelength limit (q << A, where A is a cut off) we get 

n^k), n_- q/ (k')] = -5 u #5fj?.Vno(k) + ^n (k)5._- g ,q.V5 ng . (4.3) 

We are interested in the behavior of systems with a Hilbert space restricted to the 
vicinity of the Fermi surface. We define the Hilbert space to be the filled Fermi sea \FS > 
and the tower of states obtained by acting finitely on it with local fermion operators. More 
specifically we consider a shell of states of thickness D (measured in units of momentum) 
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around the Fermi energy. Next we imagine dividing up the Fermi surface in patches |T0|JTT 



centered at points k on the Fermi surface of the filled Fermi sea. Each patch has thickness D 
and width A ( in momentum units and tangent to the surface). Instead of the sharp operators 
riq-(k) we will considered operators smeared over each patch. We will replace eq (4.3) by a 
weaker identity valid in the restricted Hilbert space. Furthermore, we will make the explicit 
assumption that we are both in the thermodynamic limit (the momenta form a continuum) 
and that the Fermi surface is macro scopically large (g < D C A < pp). 

Since the vectors k are at the Fermi surface and q are very small, in the restricted Hilbert 
space it is possible to replace the r. h. s. of (4.3) by its expectation value in the filled Fermi 
sea \FS > , namely, 

n (k) -> (n (k)) = Q(ji - eg) 
Vn (k) - V(no(fc)> = -v % 5(fi - e % ) (4.4) 

where // is the chemical potential, eg is the one-particle fermion spectrum (from which the 
Hilbert space is constructed) and vt = Veg the velocity of the excitations. The terms ignored 
in (3) vanish clS clS the size of the Fermi surface diverges. More generally, all corrections 
vanish if the limit q <C D <C A <C pp is satisfied. These are exactly the same assumptions 
that enter in the construction in one-dimension (in appendix B we prove the equivalence of 
our construction and the well known bosonization in one dimension). Notice that the state 
\FS >, which is used to normal-order the operators, is not necessarily the ground state of 
the system of interest (as in the one-dimensional case). However, this approach will succeed 
only if the true ground state belongs to the restricted space. 

Hence, within this approximation the commutators of the operators n^(k) become c- 
numbers, namely, 

n^k),n_-,(k')] = 8 Uy g8 ? -,q.v % 8{n - eg). (4.5) 

where we drop the last term in the r.h.s. of (4.3) because its matrix elements near the Fermi 
surface are down by powers of — . We now define the operator 
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a^k) = n^k)e(sgn(q)) + n^(k)Q(-sgn(q)) (4.6) 

where the sgn(q) is +1 if q points outside the Fermi surface and —1 if it points inside the 
Fermi surface at the point k. The adjoint of (4.6) is simply, 

al(k) = n^k)e(sgn(q)) + n^k)e(-sgn(q)) (4.7) 

where we used, from the definition (4.2), that nUk) = n_^(k). 

It is important to notice that, by construction, the operator defined in (4.6) annihilates 
the filled Fermi sea 

a^k) | FS) = 0. (4.8) 

Moreover, from the commutation relations (4.5) we easily obtain, 

a^k), at{k>)} =| q.v % \8{p,- e^p + 6^) , (4.9) 

all other commutators vanish. 

Eq. (4.8) and (4.9) show that the elementary excitations have bosonic character, particle- 
hole pairs, and they are created and annihilated close to the Fermi surface by these operators, 
moreover, they span the Hilbert space of low energy. But we need an interpretation for these 
excitations. It is natural now to define the coherent state JZD] 

| Uir (k)) = U{k) | FS) (4.10) 

where 



U(k) = exp £ u q ik) a\{k) (4.11) 

V q / I 1- V k I / 

where the sum (and all other sums that follow are restricted to q << A). Observe that 
from the definition (4.7) we have a_^(k) = aUk) and we choose u^ik) = u^(k). Using this 
property and the commutation relation (4.9) we find 

U-\k) a q ik) U(k) = a q ik) + 5(jji - e^v^k) (4.12) 
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which, together with (4.8), leads us to the eigenvalue equation, 



a^k) | u^k)) = 6(jj, - eAvtuJk) | u^k)). 



(4.13) 



It is easy to see that u^(k) is the displacement of the Fermi surface at the point k in 
the direction of wg. Indeed [^TJ, suppose we change the shape of the Fermi surface at 
some point k by an amount u(k). The occupation number changes up to leading order by 
5(n (k)) = — ^q^ u(k) = v^5{ji — e^)u(k) which is precisely the quantity which appears in 
(4.13). Hence, the coherent states of (4.10) represent deformed Fermi surfaces parametrized 
by the bosonic field u^(k). 

Next we define a many-body state which is a direct product of the coherent states defined 
above, 



I M) = H®U(k) | FS) = E[u] | FS) 

k 

where, due to the commutation relation at different fc's, 

/ 



(4.14) 



n[u\ = exp 



(4.15) 



The adjoint is simply, 



E,'[u] = exp 



E 

\k,q 



2 | q. Vi | 



(4.16) 



/ 



From the above equations we obtain the overlap of two of these coherent states 



({w} | {u}) = (FS | 5t[«;]S[ii] | FS) = exp £ 



k,q 



2 I q-vz 



W*Jk) Uq*(k) 



(4.17) 



It is also possible to find the resolution of the identity for this Hilbert space, 



n 

k,q 



. 2tt I q.vt | 



duJk) du^(k) | Uq~(k)) (v,cf(k) 



exp 



E 



V l 5 ^~ e k) 

2 I q.vt | 



Uq{k) 



(4.18) 
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and we conclude, as expected, that they are overcomplete p2 . 



In order to study the dynamics of these modes we can construct, from (4.17) and (4.18), 
a generating functional as a sum over the histories of the Fermi surface in terms of these 
coherent states in the form Z — J D 2 u exp{i S[u]} where S is the action whose lagrangian 
density is given by (% = 1), 

L[u] = £ V ^:jf iu 'Jk, t )^M - «»> [ H j <«» (4.19) 
11 jrt 2 | q.Vf « y ' ' dt ({«} | {«}) 1 ' 

where H is the hamiltonian written in the restricted Hilbert space which will be studied in 
the next section. 

V. THE BOSONIC HAMILTONIAN 

Consider a fermionic system described by the an hamiltonian of the form H = K + U 
where, 

K = H e ? 4°p> ( 5 - 1 ) 

V 

is the non-interacting term and 

U = — f- - ->, , - ct g c~ jet g-c-, , , (5.2) 
/'•/'''/' 

is the interaction between the fermions. We will assume that the interaction between the 

— * 

fermions is important only close to the Fermi surface. This means that the vectors p and p 1 
are at the Fermi surface in what follows. 

We see that the interaction is given in terms of the bosonic operators as, 

^=^E W' +( f n- q ip)n q ip>). (5.3) 
p,p',<f 

Of course the free part, K, can not be written directly in terms of the bosons since it 
is quadratic from the beginning. However, since we are in a restricted Hilbert space what 
really matters is the effective dynamics that it can generate. We know that the operators 



n^{k) generate the restricted space therefore it is natural to look at commutation relation 
between K and n^ik). Up to first order in q we find 

K,n^k)]=-q.vtfi3(k), (5.4) 

— * 

observe that in order to keep the vector k at the Fermi surface we have to multiply the above 
expression by a distribution function of the form 5(/i — e%)/N(0)V where 

tf(°) = 7E^- c *) ( 5 - 5 ) 

k 

is the density of states at the Fermi surface. 

Now we notice that there is another operator which gives the same commutation as in 
(5.4) in the restricted Hilbert space, namely, using (4.5) we find, 



J2 n -q'(P) n AP)' n ?(k) 



-2q.v % n^k)5{n-e % ), (5.6) 



which means that in the restricted Hilbert space we can rewrite, 

Using the definition (4.6) and the equations (5.3) and (5.7) we conclude that, in the 
restricted Hilbert space, the hamiltonian for the interacting fermions in simply given by, 

k,k',q,sgn(q)>0 

where 

S" " 

G \k' = AT(0) + f ( 5 - 9 ) 

which describes a free, quadratic, theory. 

It is worth noticing the similarity between (5.8) and (5.9) with Pomeranchuk's expression 
(2.8). This resemblance is expect since we are looking at the same kind of fluctuations. We 
would say that we have a quantum version of Pomeranchuk's construction. Moreover, all 
the scaling arguments of the section II can be used here in other to understand the stability 
of the Fermi liquid as a renormalization group fixed point. 
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VI. FERMI LIQUID PROPERTIES 



This section is devoted to show that we can get all the results of the Fermi liquid theory 
from the point of view of the bosons. In particular, the relationship between the bare mass 
of the fermions and the effective mass of the Landau theory is discussed in appendix C. We 
start by studying the thermodynamics (imaginary time) of the system and then we study 
its dynamics (real time) from the point of view of our path integral. 

Firstly we will rewrite the lagrangian in (4.19) using the hamiltonian we have just derived 
in the last section. Namely, we substitute (5.8) in (4.19) and use the eigenvalue equation 
(4.13) in order to get, 

L[u] = E Vl %~^ u ^ ( iS ^l - \ WfaSQ* - cp)) u q ik>, t) (6.1) 

k,k',q,sgn(q)>0 V ' k ^ ' ' 

where 

Observe that the vectors k and k 1 in the above expressions are restricted to the Fermi 
surface due to the Dirac's delta function. We therefore parametrize these vectors by the 
solid angles VL and Q! , respectively. We will assume throughout this paper that the density 
of states is a smooth function across the Fermi surface (which means that there is no Van 
Hove singularity present). In our language we make the following substitution, 

E ^ - H)h = J (6.3) 

k d 

where is the integral over the whole solid angle, Sd = J dQ. Furthermore, using the fact 
that the density of states is smooth we can replace the velocity by 

H = v F^k ( 6 - 4 ) 

— * 

where is an unit vector perpendicular to the Fermi surface at the point k. 
With these substitutions we can rewrite the lagrangian (6.1) as, 
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Sd j J - / 



where 



° (i<W J: - ^-(fi, SI')) u ? (Q') 



dt 



(6.5) 



F<f(fi, ft') 



(6.6) 



Observe that the field -u,j{ft) has dimensions of momentum, as required by our interpre- 
tation, however, we can rewrite the action in terms of a new dimensionless bosonic field, 
0,f(ft), by a simple scaling of the original field as, 



v F N(0)V 



1/2 



^2^1^(0)1 
and rewrite the generating functional as 



(6.7) 



*=/■■■/ n ( 

q,Cl,t,sgn(q)>0 \ 



where 



d<p q in,t)d<pt(n,t) \ iS[4>{ 
2tt r 



(6.8) 



s = Jdt(jdn Yl -Jdnjdn' Yl <W M "( Q > "0^)1 > 

V q,sgn(q)>0 ' q,sgn(q)>0 / 

(6.9) 



where 



M«(ft,ft') = J|g.n(fl)||g.7i(n0 



's d 6(n-n') 



S d V.-"v--/n.-v--/. ^ iV(0) 

We can split the field 0<f(ft) into real and imaginary part as, 
and we trace over the imaginary part in order to get, 



+ /(ft,ft') 



(6.10) 



(6.11) 



Z = {detM)~ 1 ' 2 J ... J J] 



q,Q,t,sgn(q)>0 



IT 



1/2 



dip$(n, t) 



(6.12) 



where, 
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S eff = J dtjdttj dtt' 



q,sgn(q)>0 



dt 



dt 



(6.13) 



We proceed further and expand the action in its eigenmodes by rewriting, 



(6.14) 



where 



(6.15) 



Equation (6.15) can also be rewritten in terms of the original fields if we use (6.7) and (6.10), 
v F q.n(Q) ^«J(n) + / ^F(Q,Q')4(fi')) = ^(Q). (6.16) 
Using the completeness of the states defined in (6.15) we rewrite the generating functional 



as, 



Z = (detM)-^ J ... J J] \Q 1/2 d^(t) 



where 



'^(t)\ 2 



(6.17) 



\,q,sgn(q)>0 q 



dt 



(6.18) 



which, as expected, is the action for a free massive bosonic field. 

A. Thermodynamic Properties 

The extension of the lagrangian (6.5) to imaginary time (i. e. finite temperature) is 
trivial. We will consider first the case in which only the diagonal term in eq (6.6) is present, 
namely, 



and therefore it is straightforward to show that the partition function is given by, 



(6.19) 
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Z = exp(-(3F) = || 

Q,q,sgn(q)>0 



sinh (W) 



(6.20) 



where f3 — 1/T is the inverse of temperature (our units are such that = h = 1)- 

Using the expression for the free energy, F, given in (6.20) we easily calculate the specific 
heat, 

4 n,<f, Sff n( ? )>o smh (^y^J 

We have to be careful in order to understand the integrals in (6.21). We built a reference 
frame which runs over the surface (a Fresnel frame as in differential geometry). The Fermi 
velocity (normal to the Fermi surface at pp) and gradient of the Fermi velocity (which spans 
the vectors in the plane tangent to the the Fermi surface at pp) form a local orthogonal basis 
for the vectors q. We therefore split the integral in (6.21) in a part which is normal to the 
surface, qx, and another which is in the plane tangent to the surface, qp. In the absence of 
interactions the dispersion is given by (6.6) as, 

E{q N )=v F q N . (6.22) 

Firstly we observe that the density of states is given as an integral over the tangent 
component. Indeed, from (5.5) and the definition of the local frame, we have, 

N(0) = | / - = | / <Mn)£!M. (6.23) 
b d Jv F b d J J Vp 

We also have assumed that the normal component is not affected by the curvature (see 
(6.4)). Thus we can integrate out the tangent component in (6.21) using (6.23). The final 
result is, 

— = T N(0)v F J Q ^ g , nh2(MfM) . (6.24) 



Using (6.22) we finally find, 



% = yiV(0)T (6.25) 
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which is the expected result for a Fermi liquid ||15|| . Observe that the dimensionality plays 



no role here except in the calculation of the density of states. 

We can clearly see the differences between the bosonic excitations described in this paper, 
which are due to coherent superposition of particle-hole pairs, and the usual free bosonic 
excitations. First of all, naively we would expect that the specific heat should behave as 
T 3 as in usual free bosonic theories since the bosonic hamiltonian (5.8) is quadratic in the 
bosons. However, the the bosonic field in the Fermi liquid lives on the Fermi surface, that 
is, it is topologically constrained. The Fermi velocity , through the density of states, defines 
the metric of the manifold where the fields propagate. In the case without interactions the 
bosonic fields oscillate without coherence (like decoupled harmonic oscillators with phases 
distributed at random). These oscillation are responsible for the contribution of the particle- 
hole continuum to the specific heat as shown in (6.25). 

We can proceed further and evaluate the effect of the interactions in this problem. The 
general approach would be to calculate the eigenvalues of equation (6.15) and evaluate the 
bosonic determinant in the partition function. However, we use an approach which is inspired 
by the calculations of the Fermi liquid theory Firstly, we expand the interaction term 
for small angles by defining the vector p = k' — k such that 

h,w = a + b (^.^j (6.26) 

where a and b are parameters which depend on the specific form of the interaction. Due to 



the geometry of the interaction we can rewrite (6.26) as [11 



fan, = o + 6 2 f (6.27) 

Since we are working with a local expansion we will replace the non-local term in (6.6) 
by a local one which is an average over the surface, 

/n,n' -> ( I dn'f n ,n>) <W -> 2 f N 2/ M n>, (6-28) 
V J y q z N + p 2 /4: 

changing the sum to an integral and taking into account that this integral must be done on 

the surface we found that the change in the energy due to interactions can be written as 
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AE(q N ) = / dp-f—^ (6.29) 

where A is the cut-off. In three dimensions, we have, 

8bv F N(0) 3 fq N \ 

E(q N ) = v F q N 2 W n ~T ( 6 - 30 ) 

P F V A / 

and in two dimensions one finds, 



E(q N ) = v F q N + nb^^lq 2 N . (6.31) 

Pf 

Observe that these expression depend only on the normal component and we can use 
equation (6.24) in order to calculate the corrections for the specific heat due to interactions. 
Up to first order in b we have, 

6Cv -aTHn ( (6.32) 



V \Av F 



in three dimensions, and, 



SC V _ 2 



7]T 2 (6.33) 



in two dimensions, where, 



a = 15^4 (6 ' 34) 



and 



V= 3 — ' 6 - 35 

71V F 

where ((n) is the Riemman Zeta function of n. 

The result (6.32) is well known |[T5|p3| (apart from the prefactor which is different due 



to the average over the Fermi surface we have carried out in (6.28)) and the result (6.33) 
was obtained recently using the RPA approximation |24 . 



We conclude therefore that in our approximation we are counting the correct number of 
states at low temperatures. This result proves the consistency of our method. 
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B. The semiclassical dynamics 



Notice that the lagrangian (6.5) is quadratic in the fields and, therefore, the semiclassical 
approximation is exact. The semiclassical equations of motion for these lagrangian are given 
by the saddle point equation (the Euler-Lagrange equations) derived from L, namely, 



i- 



ot 



— F{Vt,n')u^t',t) (6.36) 
Sd 



where we have defined the angle 9 by q.n(Q) = cos 9. Observe that if, 

u g in,t) = e-^u^tt), (6.37) 

we recover equation (6.16). 

Eq. (6.36) is the Landau equation of motion for sound waves (the collective modes) of 
a neutral Fermi liquid where ft g is the scattering amplitude for particle-hole pairs |jT3| . 
Observe that the Landau equation gives the eigenmodes of the problem (see (6.15)). 

We conclude therefore that our bosons are the sound waves which propagate around the 
Fermi surface at zero temperature. The solution of (6.16) (or (6.36)) will give the possible 
values for the frequencies of oscillation for these modes and they will depend essentially 
on the Landau parameters of the theory. The Landau equation, (6.16), yields solutions 
which represents both stable collective modes ( "sounds" ) as well as solutions with imaginary 
frequency which represent the particle-hole continuum. This behavior is a direct consequence 
of the phase space. Notice that in one dimension these unstable solutions are absent and 
only the collective mode is left. 

In order to illustrate the behavior of these sound modes we present here an explicit 
calculation in two dimensions. This calculation will enable us to understand many interesting 
features of these sound modes. 

As in the section II, we expand the interaction and the displacement in Fourier compo- 
nents, 

oo 

F(6, 9') = F n e l(9 - e ' )n (6.38) 

n=— oo 
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where F_ n = F* = F n due to the symmetry of the problem and, 

oo 

<0)= E u ne Wn , (6.39) 

n=— co 

where w_ n = u* n . 

From (6.16) we have, 

r 2n dO' 

(s - cos 6)u(9) = cos 9 —F(6,e')u(6') (6.40) 

Jo 2tt 

where vps = u/q is the velocity of the sound waves. 

Using (6.38), (6.39) and the orthogonality relation of the Fourier components it is easy 
to show that (6.40) can be written as a matrix equation, 

CO 

(l + F n )u n = ]T F m K(m-n)u m (6.41) 

m=— oo 

where 

K (n) = f 277 ^. s -—e me . (6.42) 

v ; Jo 2ns-cos9 v ; 

We allow s to have an infinitesimal imaginary part and we integrate (6.42) in the complex 
plane in a contour in the unit circle around the origin. We found 

K(n) = -^L=( sgn (s)r (s - N \s\ > 1 

= . S (s-Wl- s 2 ) H \s\<l. (6.43) 

iy/1 — S 2 v ' 

It is easy to see from (6.43) that all the modes are damped if the velocity is smaller of 
the Fermi velocity, \s\ < 1. This result is clearly expected, it is nothing but the Landau 
damping of the collective modes [plj . 



Here we are only interested in the stable solutions of (6.41). Therefore, for s > 1 we can 
replace (6.41) by 

CO 

(1 + F n )u n = F m K(m-n)u m (6.44) 

m=— co 

where we define a kernel, 
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K(n) = A(s)e- lnla{s) (6.45) 

with, 

A(s) = (6.46) 

and 

a(s) = In (s + Vs 2 - l) . (6.47) 

For a purely local interaction in the real space the interaction is angle independent (long 
range over the surface in momentum space) and we define, 

F n = F S nfi . (6.48) 

Substituting (6.48) in (6.44) for the n = 0we obtain the allowed value of s, 



-TO 



'('Hi -' 

and for n^Owe obtain the value of the Fourier components 



(6.49) 



u n = F K(n)u . (6.50) 

Substituting (6.50) in (6.39) and summing up an ordinary geometric series with a little help 
of (6.49) we easily get, 

u(6) = u F (6.51) 

s — COS V 

which is expect result if we have just solved (6.40) for a constant interaction. The result 
(6.51) is the well known zero sound mode which has the same shape as the zero sound in 
the three dimensional case [ |I5|j . 

We can explore even more the matrix equation (6.44) due to the simplicity of the kernel 
in (6.45). We are mainly interested in the behavior of the Fermi surface in the forward 
scattering direction since it has been argued || that it is in this direction that pathological 
effects can happen. 
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Since we are not interested in the behavior over the whole Fermi surface we will transform 
the matrix equation (6.44) in a simple second order differential equation by expanding the 
kernel in terms of the Fermi wavelength. 

Consider the arc-length defined on the Fermi surface by, 

v = p F 9, (6.52) 

the expansion (6.39) is rewritten as, 

u{y) =Y,e iXu u{\) (6.53) 

A 

and due to the periodic boundary conditions, u{y + 2irp F ) = u(v), we must have, 

\ n = — . (6.54) 
Pf 

Observe from (6.54) that in the limit of p F — > oo the matrix equation (6.44) can be 
replace by an integral equation, 

/+oo 
d\'G(X - X')F(X')u(X') (6.55) 

where 

G(X)=2- — — 2 (6.56) 

a J-oo Z7T y _|_ [ x ) 
\PFaJ 

where we have used (6.45) and an integral representation for the exponential function. Ex- 
panding (6.56) up to second order in p F l one finds, 

If we substitute (6.57) in the original equation (6.55) we find a differential equation which 
should be valid for small arc- lengths compared to 2-np F . We proceed further and define the 
following quantity, 

z(X) = F(X)u(X) (6.58) 
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which, from the arguments given above, obeys the differential equation, 

+ U{s,\)\ z{\) =e{s)z{\) (6.59) 



where 



= — F\X) — ^ ' 



and 



= (pMs)) 2 \2^& - lj . (6.61) 

Eq. (6.59) has the form of a time independent Schrodinger equation for a particle in a 
potential U with energy e. Once we know the boundary conditions for the problem we can 
determine its eigenf unctions (and therefore we determine u using (6.58)) and its eigenenergies 
(and we determine s using (6.61)). 

In this paper we will assume a particularly simple form for the interaction term which 
has been used in the literature in order to investigate the problem of interactions of finite 



range p5" | 



F{y) = 7r 7 p F A c e- AcM . (6.62) 

There are two parameters here, 7 gives the strength of the interaction and A c plays the 
role of the range of the interaction. When A c —>■ the interaction is angle independent 
(is short range in the real space) and we obtain the zero sound solution already discussed; 
when A c — > 00 the interaction is long range in real space and strongly angle dependent in 
momentum space. This parameter will enable us to study the crossover from short range to 
long range interaction in terms of the dynamics of the sound waves. 
The Fourier transform of (6.62) is given by, 

*W = 75^55, («■«) 

which is a lorentzian. 

30 



Using (6.63) in (6.59) we obtain the equation for the one-dimensional harmonic oscillator, 

/ d 2 



d\ 2 



+ (3%s)\ 2 z(X) = E(s)z(X) (6.64) 



where 

1/2 



and 

,1/2 



E(s) = e(s) - I ^ ) . (G.66) 



Of course the natural boundary condition is that the function z vanishes at the infinity. 
The solution of (6.64) is standard. The eigenf unctions are, 

z n {\) = e-^ x2 H n {(3 2 \) (6.67) 

where H n is the Hermite polynomial of order n and the eigenenergies are, 

E(s n ) =2(3 2 (s n ) (n+l). (6.68) 



2, 

Equation (6.68) is self-consistent and it gives the allowed values of s. If we substitute 
(6.46), (6.47) in (6.68) we get the following transcendental equation, 



l + -)a(s n ) + 2 
7 



1 ' n +\)- A (Sr, 



0. (6.69) 



J 1/2 Pf\c V 2, 

It is very easy to check that the value of s n decreases monotonically with n and, in particular, 
Urrin-tooSn = 1. This result means that the modes approach asymptotically the particle-hole 
continuum. 

Furthermore, substituting (6.67) and (6.63) in (6.58) we find 

u n (\) = (A 2 + X 2 c )e^ 4x2 H n ((3 2 X). (6.70) 
We now Fourier transform (6.70) back to the arc-length in order to get, 

«.(,)= (2„ + l + ^„)A c -(^) 2 )e-«^)V„(^). (6.7!) 
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Equations (6.69) and (6.71) are the solutions of (6.44) in the limit when the Fermi 
momentum is much larger than the momentum transfer between the particle-hole pairs 
(which means, as we said, that we are looking close to the forward scattering direction). 

We want to change the range of the interaction, A c , but we keep the density of the system 
constant. In this case the product PfK is finite ( X c Pf ~ ^ where r is the mean distance 
between the particles and is the Bohr radius). In other words, it means that (5\ c is finite 
in what follows. Therefore, keeping the above product constant and sending A c — we find 
from (6.65), 

(3 — oo (6.72) 

and from (6.71), 

u n {v) — constant, (6.73) 

if n is even and u n {y) = if n is odd. Observe that this result is in agreement with the zero 
sound result, (6.51), which is constant in the forward direction (8 ~ 0). 
However, for A c — > oo it is easy to see that, 

/3 - 0, (6.74) 

and 

uM ~ (i) c ~* ( * )aj3 » [fj ^ 

which oscillates strongly near v = 0. We would argue therefore that, while for short range 
interactions (A c — > 0) the field u(u) is smooth around the forward scattering direction, for 
long range interactions (A c — oo) the field becomes rough near the same direction. This 
could represent a signal that there is something particularly non-conventional happening 
there. Notice, on one hand, that this result has similarities with the problem of long-range 
order in low dimensional systems |2Sj since it is well known that one- dimensional systems 



at zero temperature with local interactions do not have long range order. However, on the 
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order hand, also quantum one-dimensional interfaces at zero temperature can not have long 



range order while two-dimensional interfaces always have long range order nj. Since the 
Fermi surface has the properties of a real interface in momentum space it could be argued 
that the Fermi surface could undergo a roughening transition of the Kosterlitz-Thouless 
type for singular interactions. 

However, the calculations we have carried out are purely classical and up to now we do 
not know if quantum fluctuations can wash out this behavior. It is worth mention that in 



quantum crystals the quantum fluctuations can wash out the classical ones JT4 



C. Quantum fluctuations of the Fermi surface 



It was shown some time ago |26| that quantum one-dimensional crystals can lose long 
range order when the interactions are sufficiently local. Quantitatively it means that the 
equal-time correlation function for the deviations of the equilibrium position of the atoms, 
u(R), diverges logarithmically for long distances, that is, (u(R)u(R')) ~ ln\R — R'\ as 
\R — R'\ — > oo. We would argue, due to the classical calculations we have carried out 
that the same could happen for the fields Ug-(Q) for long range interactions in real space. 
Therefore, in order to understand this problem we calculate the quantum correlation function 
for different pieces of the Fermi surface. 

From (6.18) is straightforward to calculate the correlation function for different pieces of 
the Fermi surface, namely, 

(T^{fi,t)^(fi',t')> = E^K'( fi ') (T^(t)^(t')) (6.76) 

A, A' 

where T is the chronological operator and v^(Q) are the solutions of (6.15). The correla- 
tion function in the r.h.s. of (6.76) is easily calculated from (6.18) (we assume that the 
eigenfrequencies are well behaved functions of q and A) and it reads, 

(T^mp')) = V'^ e ^ M '' (6.77) 
Thus from (6.76) and (6.77), 
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(r^n^o^O) = <W E^K^O e ^ M '' ( 6 - 78 ) 

A 



and therefore the equal-time correlation function is simply (use the completeness of the 
states in (6.15)), 

(T^(n,t)^(n',t)) = 6^6{Q - n'). (6.79) 

Eq. (6.79) shows that the correlation function is always finite and therefore there is no 
roughening (there is long range order), that is, the Fermi surface is always smooth. 

The comparison between this result and the result of quantum crystals |T3J we know that 



this result is due to a conservation law, namely, the local conservation of the volume of the 
Fermi surface (the Luttinger's theorem ||13|| ). From the other point of view this is expect 



since the sound waves obey a diffusion like equation, (6.36), and therefore a bump on the 
Fermi surface diffuses leaving nothing behind it. 



VII. CONCLUSIONS 

We construct a bosonization of a Fermi liquid in any number of dimensions in the limit 
of long wavelengths. The bosonization is valid for a restricted set of states close to the 
Fermi surface. We generate a set of creation and annihilation operators which span all the 
restricted Hilbert space. We showed that is possible to construct a set of coherent states 
(bosonic shape fields) which are coherent superposition of particle-hole excitation close to 
the Fermi surface and which are the displacements of the Fermi surface in some direction. 
The physical interpretation for the existence of these fields comes from the observation that 
the Fermi surface is a real quantum object which is responsible for the whole physics of Fermi 
liquid systems at low energies. From the coherent states we generate all the thermodynamics 
and semi-classical dynamics of a Fermi liquid. 

We have shown that an interacting hamiltonian of the original fermions is simply 
quadratic in terms of the bosons, which means that the bosons are free and move on the 
Fermi surface like sound waves. They are not free fields in the usual sense since they are 
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topologically constrained to live in the Fermi surface. This new feature produces very inter- 
esting results. And, in particular, reproduces all the results of the Fermi liquid theory from 
a bosonic point of view. 

We obtain the correct thermodynamics of a Fermi liquid even with the corrections due to 
the scattering between the fermions. This observation lead us to conclude that our approach 
is consistent and that we are counting the correct number of states in the Hilbert space. 

Moreover, we have shown that the semiclassical dynamics of these sound waves is de- 
scribed by the Landau theory. In two dimensions we calculate the form of these sound waves 
explicitly for short range interactions (zero sound) and long range interactions. We have 
shown that, at classical level, we could conclude that it is possible to have a roughening 
like transition of the Kosterlitz-Thouless type, which is closely related with the absence of 
long range order in one dimensional crystals at zero temperature. However, we show that 
quantum fluctuations destroy this effect due to the local conservation of the volume of the 
Fermi surface (Luttinger's theorem). 
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APPENDIX A: FERMI LIQUID AND NON-FERMI LIQUID BEHAVIOR 



In this appendix we discuss the possibility of other behaviors for the imaginary part of 
the self energy than the conventional Fermi liquid one, (3.13). We will assume yet that the 
real part of the self energy is linear in the frequency close to the Fermi surface but for the 
imaginary part we will assume that, 

Ej(jfe, w) = C % \w- sgn(fx - w). (Al) 

If we substitute (Al) in (3.12) and make the same change of variables which lead to 
(3.14) we find, 

E S = /i-iiv^i^ (» + el--\5k\\V6 f < ----- ^ 



^ 2tt k x PF 7 (zfx 2 "{x + sgn(5k)) 2 + C\5k 2 " | Ve ?F \ 2v )' 

(A2) 

When v > it is easy to calculate the limit analogous to (3.15), 

^™^o y— — - *' - ' ^ ' - — = tt | Z & | 5 (*"(* + s^)) , (A3) 

[ZJx 2v [x + sgn(Sk)) 2 + C\5k 2u | Ve^ | 2 ^J 

which gives the same result as in (3.16) (Fermi liquid behavior). However, if v < the above 
limit takes a different form, namely, 

^° (z:V"(x + S ^n(^)) 2 + C\5k 2 » | Ve p> | 2 ) C % \8k\" \Ve%\^ 



and therefore, substituting in (A2) we find, 

e s , = ( - 1 - cr 2 \6k\v4 \\ {u ~ 1)2 ) e(-Sk). 

The discontinuity is easily calculated in the limit as 5k — > 0, 



(A5) 



AS(pf) = I P> 2(1 !^ } I l^l ( ^ 1)2 - 0, (A6) 

and therefore there is no singularity left (non- Fermi liquid behavior). 
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APPENDIX B: BOSONIZATION IN ONE SPATIAL DIMENSION 



In this appendix we show that the formulae obtained in section IV of our article are 
consistent with the bosonization procedure in one dimension. In order to do so we prove 
that we can obtain the same commutation relations for the densities as it is obtained in the 
standard procedures of bosonization. Therefore, all other results relative to the calculation 
of the fermionic operators and the one particle Green's function follows from this result. 

We will concentrate particularly on the Tomonaga's model PJ but the extension of the 
arguments to the Luttinger's model HQ is absolutely straightforward. 

The Tomonaga's model for spinless electrons is described by the following hamiltonian, 

# = E ^ v c v + V q p(q)p(-q) (Bl) 

where the dispersion relation is given by, 

e p = v F \p\, (B2) 

V q is Fourier transform of the electron-electron interaction, L is the length of the system 
and 

P(9) = E%(*0 (B3) 
fc- 
is the density operator. 

Now we split the density operator in two terms, one for right movers and other for left 
movers, 

k>0 

p 2 (g) = £n 9 (fc), (B4) 

fc<0 

thus, p(q) = Pi(q)+p 2 (q)- 

The one dimensional version of the commutation relation (4.5) is, 

[riq(k),n- q >(k')] = sgn(k) 5 k #5 m < qv F 5(p - v F \ k |). (B5) 
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Therefore we have from (B4) we easily find, 



[p l (q),p 1 (-q')] = (-l) J+ \j<lVf E S(n-v F \k\). 

fe(-l)J<0 

where the chemical potential is written as fi = vfPf- 

In the thermodynamic limit we transform the sum into an integral (J2k 
we finally get, 

-iQL 



2tt 



i,3 = l,2 



(B6) 



^ j dk) and 



(B7) 



as expected for spinless fermions JE]. 

The bosonic operators are defined in terms of the densities for q > as (compare with 
(6.7)) 

I2n 



with their adjoint they obey the usual bosonic algebra, 



(B8) 



b k , b\, 



(B9) 



Following the same kind of argument of section V we can prove that the Tomonaga's 
model is purely quadratic in terms of the bosons and can be easily diagonalized by a Bo- 
goliubov transformation. In particular the calculation of the correlation functions follows 



exactly as in the work of Mattis and Lieb [|J and Luther and Peschel for the Luttinger's 
model. 



APPENDIX C: THE EFFECTIVE MASS 

In this appendix we show the complete consistency of our approach and the Landau 
theory of the Fermi liquid via the calculation of the effective mass of the Landau theory in 
terms of the bare electronic mass. 
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Suppose we displace all pieces of the Fermi surface in momentum space by the same 
infinitesimal amount, q. This is equivalent to look to the system from a referential which 
moves relative to it with constant velocity. By simple geometric arguments it is easy to see 
that the displacements of the Fermi surface will change by the quantity, 

u q ik) -> u g ik) + (CI) 
v k 

The total change in the action can be obtained from (6.1). Up to first order in q one 
gets, 



k,k',q '' h 



Using (6.2) and (5.9) we can rewrite the above expression as, 



^ = - E WjE^O* - > (C3) 



where, 



Notice that since we are transporting the whole Fermi sea by a constant vector the change 
in the lagrangian must be minus the total change in the energy of the system. Therefore, 



is the contribution for current carried by a fermion with momentum k |L5| . Moreover, if the 
system is homogenous and isotropic the current and the velocity must point in the direction 
of k. Naturally, the current transports a mass equals to the mass of the fermion, m, namely, 

h = -• (C5) 
m 

In the Landau theory the effective mass, m*, is defined as the coefficient of the velocity, 

v % = A. (C6) 
m 

Substituting (C5) and (C6) in (C4) and assuming k and k' at the Fermi surface and 
k.k' = pp cos 9 we easily get, 
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777* 1 « „ 

— = l+yE/** COB* (C7) 

k' 

which is the usual relation between the effective mass and the bare mass in the Landau 
theory of the Fermi liquid [|1^ . 



In three dimensions we can expand the interactions in Legendre polynomials, 

F(#)=^F L P L (cos#) (C8) 

L 

where F = N(0)f. Use the density of states as defined in (5.5) and the orthogonality 
between the Legendre polynomials in order to get, 

- = 1 + ^. (09) 
m 3 

In two dimensions we use (2.6) and the above definitions and we find, 

^ = 1 + *. (CIO) 

The galilean invariance present in the argument of this appendix shows that the mass 
which appears in the Fermi velocity in the case of a homogenous system is not the bare mass 
but the effective mass which is calculated in the above expressions. 
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